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Abstract. The Sunyaev-Zel'dovich (SZ) effect has a distinct spectral signature that 
allows its separation from fluctuations in the cosmic microwave background (CMB) and 
foregrounds. Using CMB anisotropics measured in Wilkinson Microwave Anisotropy 
Probe's five-year maps, we constrain the SZ fluctuations at large, degree angular 
scales corresponding to multipoles in the range from 10 to 400. We provide upper 
bounds on SZ fluctuations at multipoles greater than 50, and find evidence for a 
hemispherically asymmetric signal at ten degrees angular scales. The amplitude of 
the detected signal cannot be easily explained with the allowed number density and 
temperature of electrons in the Galactic halo. We have failed to explain the excess 
signal as a residual from known Galactic foregrounds or instrumental uncertainties 
such as l//-noise. 
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1. Introduction 

The Sunyaev-Zel'dovich (SZ) effect [l| is the inverse Compton scattering of the cosmic 
microwave background (CMB) by electrons throughout the universe. The SZ effect can 
be partitioned into two components: the kinetic effect (also known as the Ostriker- 
Vishniac effect ^) due to bulk motion of electrons with respect to the rest frame of 
the microwave background, and the thermal effect due to energy transfer from hot 
electrons in massive galaxy clusters The SZ thermal signal is expected 

to be important in future constraints on the underlying cosmology of the universe, by 
revealing the properties of clusters in cluster-counting experiments (e.g. see Carlstrom 
et al. 2002 jol) and via the angular power spectrum 10 . 

At degree angular scales, there are no detections of the SZ fluctuations as primordial 
anisotropies of the CMB dominate the SZ signal by at least 3 orders of magnitude. 
Fortunately, a separation of the SZ anisotropies from those of the CMB is possible 
due to the fact that the SZ effect has a distinct frequency spectrum, as the inverse- 
Compton scattering on average increases the net energy of the CMB photons and move 
photons from the low frequency Ray leigh- Jeans (RJ) tail to higher frequencies In 
multi-frequency CMB experiments spanning a wide range of frequency coverage across 
the SZ null frequency at 217 GHz, one can separate the CMB down to sub-percent 



level required to study the SZ anisotropy power spectrum [ll|. This approach was first 
applied to constrain the sub-degree SZ fluctuations in BOOMERanG 2003 data 121 . 
Here, we analyze the WMAP five-year CMB data to constrain SZ anisotropies at degree 
angular scales corresponding to multipoles £ = 10 — 400 on the sky. We ignore £ < 10 
as a measure to avoid biases from previously reported anomalies, the divergence of 
1 / /-noise, and large uncertainties associated with simulating the CMB sky properly to 
match WMAP data at these multipoles Q, Q [isl, Q, [13, 0, H, 5 . 

In Section 2 we detail our approach to extract the SZ signal from the multifrequency 
WMAP5 data set; in Section 3 we present our results; and in Section 4 we provide a 
discussion of these results. 



2. CMB, Noise, and Foreground Removal 



We distinguish the SZ signal from other sources of anisotropies by minimizin g th e 
covariance relative to the SZ frequency dependence (also see Veneziani et al 2009 [12|). 
Although this tech niq ue is originally employed for removing foregrounds from CMB 



anisotropies 2l|, |22, l23[ , we utilize it to recover SZ fiuctuations from which the primordial 
CMB is subtracted along with other sources of noise. 

We measure the raw SZ_ power spectrum from a weighted mean of the spectra in 
different frequency bands 
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C, 



rawSZ 



wJC^w^, 



(1) 



where the scale dependent weights at each frequency v are obtained by minimizing 
the contribution of the CMB, noise, and foregrounds in the variance CY™^"^ (of multipole 
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moments oF^^"^ = ^j^^a£rn{t^i) as elucidated below). We enforce the constraint that 
SZ estimation is unbiased, via wi{vi) = 1- The optimal weights for reconstructing 
the SZ are thereby obtained from applying the method of Lagrange multipliers to Eqn. [1] 

where e is a unit vector, such that e{vi) = 1 Vz/j. 

We construct the positive-definite covariance matrix for frequency bands i and j as 
a function of angular scale, wherein the SZ frequency dependence is removed to prevent 
its minimization: 

( fi\ {^\m^\m) / /sky /q\ 

where /sky is the sky fraction analyzed, the multipole moments aim are extracted from 
either data or simulated maps at separate frequencies, and s{v) = 2 — (a;/2)coth(x/2), 
with X = hu/kTcMB ~ z^/56.8 GHz, is the frequency dependence of the SZ at the 
center of each of the WMAP frequency bands relative to the CMB. In the RJ limit 
s{iy) — 7- 1, such that Cf^i^u, z/') = s(z/)s(z/')C|^, where Cf^ is the SZ power spectrum in 
the RJ limit. For the two preferred WMAP5 bands [V,W] = [60.8,93.5] GHz, we find 
s(z^) = [0.906,0.784]. For the data, the multipole moments are obtained from WMAP 
five-year maps 

data _ WMAP /l fA\ 

normalized by the measured beam window function bfy in the respective frequency band. 
To compute the covariance matrix of non-SZ sources, we Monte-Carlo simulated 



a =et of 250 Gaussian sky maps with HEALPifl based on the CAMB-generated |2| 

input power spectrum in accordance with the WMAP5 best-fit ACDM cosmology 
([fie/^-^fi^,/i^^^A,n„r, A|] = [0.1099, 0.02273, 0.742, 0.963, 0.087, 2.41 x 10"^]). We find 
consistent results with simulations using the exact WMAP5 measured CMB power 
spectrum in each of the bands. From these frequency independent Gaussian maps, 
masked along the Galactic plane via KQ7^ (admitting 71.6% of the sky), we then 
extracted the moments up to £ = 400 with HEALPix. To account for the masking 



in our estimated spectra, we used the algorithm described by Hivon et al (2002) [25 . 

In addition to these CMB simulations, we create 250 noise maps for each frequency 
band, given by 

iV.(n) = -^n(n), (5) 

where Ni,{h) is the noise map, n(n) is a white noise map, the frequency dependent A'obs 
is the number of observations per pixel, and the noise per observation ctq = [3.133, 6.538] 
mK for the [V, W] bands§. We moreover explore the impact of l//-noise in our 
simulations. To this end, we assume similar l//-noise properties for five-year data as 

I http://healpix.jpl.nasa.gov 

§ http: / /lambda. gsfc.nasa.gov / product / map / dr3 / 
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for one- year data consistent with the observations (see Eqn. 8 and Table 1 of Hinshaw 



et ah 2003 26|), and take the same amphtude ratio between one- year data and five- 



year data of 1//- noise as the corresponding amphtude ratio of the white noise [27[. On 
the largest angular scales, the l//-noise boost in the V-band is a factor of two, and 
in the W-band the boost is at most a factor of three (see Fig. [1]) relative to the white 
noise amplitude. We comment on our result with and without 1 //-noise included in 
the discussion. From the set of noise-generated maps, we produce the corresponding 



multipole moments. 



''lm\ui 



out to the same angular scale as their CMB counterparts. 



We correct for the beam according to 



''im 



The WMAP5 maps represent a weighted mean of the five single-year maps 28| . In 



subtracting from the unreduced maps synchrotron and free-free emission, primarily at 
frequencies below 60 GHz, and thermal dust emission which dominates at frequencies 



above 60 GHz, the WMAP team 28|, |29|, ISOj utilized a temperature difference map 



of the K (22.8 GHz) and Ka (33.0 GHz) channels, an extinction-corrected H« map of 



Finkbeiner (2003) 3]J], and dust emission Model 8 of Finkbeiner et al (1999) 32 



We account for potential systematic uncertainties in the WMAP team's foreground 
removal by adding to the multipole moments of the CMB and noise simulations a 
contribution from foreground residuals. We further consider 1% uncertainties in the 



beam transfer functions of each frequency band 33|. Consequently, the multipole 



moments from simulations take on the final combined form 



''im ' "'Imu 



where B,{IL) = ht,u (l + ^h-^li^)), = 4 



sync 



;i + Arf,7^), = U{1 + /^U^l), and = 
Sy{l + Asi,7^), such that 7 is a frequency dependent, but multipole independent unless 
otherwise specified, Gaussian random number drawn for each of the 250 simulations. 
The frequency dependent dy, f^, and are the WMAP foreground template amplitudes 
of the considered dust, free-free, and synchotron emissions, respectively. These 
amplitudes are given in Table 2 of Gold et al (2009) 28|. The dimensionless Ab, 



Ad^, Af^, and Asi, respectively encapsulate a 1% beam uncertainty, and a [4.4, 



5.6]% uncertainty on the amplitudes of the dust, free-free, and synchrotron foreground 
templates in the [V, W] bands. 

In practice, we only simulate the deviation from the WMAP foregrounds and utilize 
the foreground cleaned maps in our analysis. The synchrotron foreground template was 
constructed by the difference of the WMAP5 temperature map in the K band with that 
of the Ka band, whereas the other two templates were obtained from the LAMBDA§ 
website. The ~ 5% foreground template uncertainties were obtained from the stated 

More specifically. 



28 



15 /iK uncertainty on the removal of the combined foregrounds 
we add all of the foreground templates together (multiplied by the respective amphtudes 
in each band), calculate the dispersion at of the combined template pixels, and let the 
fractional error on each foreground amplitude equal ^^J^^ ■ Table [Ushows the foregrounds 
render a bias and increase the uncertainty in the SZ estimation primarily on large scales 
{£ < 50), while the beam uncertainty is manifested on the smallest scales {£ > 150). 
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Figure 1. Top: The average simulated auto spectrum (Eqn. [3]) in the W frequency 
band for the combination of Gaussian CMB and noise illustrated by a black solid 
line about its la yellow error band. The red dotted line stems from pure CMB, and 
the orange long-dashed line is the noise contribution (beam normalized), while the 
dashed cyan line solely encapsulates the white noise. The purple and green dotted 
lines are the spectra obtained from the foreground reduced and unreduced WMAP5 
maps, respectively. The dot-dashed blue line represents the signal from the theoretical 
Komatsu-Seljak SZ power spectrum [lol |. Bottom: Optimal weights constructed from 
the covariance matrix following Eqn. [2l The upper set of dotted lines are obtained from 
the V-band and the lower set of solid lines from the W-band. The mean simulated 
weights run through the yellow bands, representing the la uncertainty about the mean. 
The foreground cleaned and uncleaned weights from WMAP5 are also illustrated for 
each frequency band. 
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In addition to our account of residual foregrounds in the error analysis, we add a 
flat spectrum of radio point sources to the temperature power spectrum of the above 
simulations, such that the numerator of Eqn. [3] becomes C|'™ — )■ _j_ ^^ps^ ^-^^ 



the frequency dependent model for the point source power spectrum 34 



where Vi and Vj are the central frequencies of the two considered bands, normalized by 
that of the Q band vq = 40.7 GHz. The point source amplitude Aps = (11.1 ± 4.1) x 
10"'^ /iK^sr for the band combination VW, and the spectral index of the point source 



flux a = -0.09 ± 0.18 34 . Moreover 



(e"" - 1)2 

K^) = (8) 



where as before x = hu/kTcuB- Beyond point source removal, we account for the 
uncertainty in the amplitude and spectral index by independent Gaussian realizations 
around the mean. In Table [H we see that point sources primarily bias our results on 
small scales. We have ignored clustering of point sources, which may add some further 
power on these small scales. 

The multipole ranges for our bins are listed in Table [H We obtain the binned 
spectra by weighting each Ce (averaged over the 250 simulations for the simulated case) 
with the respective uncertainty cr^ at that multipole. 

3. Results 

We show the contributions to the power spectrum from the CMB, noise, and their 
combination, along with the WMAP five-year spectrum and the theoretical Komatsu- 



Seljak SZ power spectrum lOj] for the W band in Fig. [T] (top). The WMAP data trace 
the simulated spectra closely, which become shot noise dominated for i > 350. 

In Table [T|, we state the binned weights of the fore- ground reduced WMAP5 data 
for the two frequency bands. The optimal weights at each multipole are illustrated in 
Fig. [1] (bottom) for both the data and simulations. The two bands have weights roughly 
symmetrically aligned along zero, and decrease towards smaller angular scales. The 
residuals are the average spectra measured on our SZ-free simulations and represent our 
bias. The total residual is slightly different from the sum of the partial residuals due to 
correlation between components. The uncertainties in the table are the 2a dispersion 
cTbin measured with our simulations. The final SZ spectrum values are corrected from 
the residual bias. The total error budget is given by the rms of the partial uncertainties. 
While in the second and third bins we find no detectable signal, in the first bin with 
10<i< 50, we find an excess signal of 58.6 ± 17.3 //K^ (95% c.L). 

For comparison, if we include a l//-noise contribution in our simulations by a large- 
scale boost to the white noise (by a factor of 2 to 3), the first bin decreases to 40.3 ±20.1 
fiK'^ (95% c.L). Removing the channels most susceptible to 1//- noise (W2 and W4 via 
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Hinshaw et al. 2003 26[) does not statistically modify the signal from that originating 
from the analysis without 1//- noise included. Unless 1//- noise is underestimated by an 
order of magnitude, it seems unlikely that excess noise can explain away the detected 
excess signal. Our findings are robust to an analysis of single-year maps in separate 
differencing assemblies (whereby Eqn. [3] is extended to an average over all combinations 
of differencing assemblies in bands i and j for all auto and cross- years) , whereas a 
cross-correlation only analysis leads to a non positive-definite covariance. 

In Fig. |2] we show our final estimate of the SZ-like power spectrum in distinct 
multipole bins for a combination of the [V, W] bands. There seems to be a clear 
non-zero signal at the largest angular scales (10 < £ < 50). Although foregrounds 
have already been reduced, this signal could potentially be due to further systematic 
uncertainties in the foreground removal. The limits on the SZ signal from a 2003 



BOOMERanG flight |12l . |35| have been included for comparison. In addition, if we 
modify Eqn. [1] from CY™^'^ = wJC^w^ to include a noise-subtraction, such that 
^rawsz _ _ ^^^^^ then our result in the first multipole bin is altered from 

58.6 ± 17.3 /iK2 to a value of 69.1 ± 18.5 /iK^ (95% c.L). 

As a simple consistency check on our calculation, we computed the leftover signal 
from the difference of the V and W maps, i.e. the power spectrum of a map given by 
v(n)-w(n) normalized by the beam (the beam is the same to sub-percent level for 
the two frequency bands at £ < 50), rendering a value of 697 /xK^ in the first bin. We 
subtracted from this quantity the mean simulated noise on the difference map, given 
by (^""^■''^^^■^^^^J^j^^^'^^''^'^^ ) = 53g ^K^. The signal derived from this method is then 
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SZ Power Spectrum Estimates 





Bin 1 


Bin 2 


Bin 3 


£-range 


10 - 50 


50 - 150 


150 - 400 


Optimal weights 


■J^eO-SGHz 


5.95 


3.69 


1.08 


"W^QS.SGHz 


-4.95 


-2.69 


-0.0831 


Raw SZ 


259 


1740 


4850 


Residuals 


CMB & Noise 


139 


1650 


4690 


Foregrounds 


57.4 


26.7 


35.5 


Beam 


-0.0414 


-0.754 


7.97 


Point sources 


3.08 


21.9 


77.5 


Total residual 


201 


1700 


4810 


SZ Band Power Uncertainties 


CMB & Noise 


15.5 


50.3 


49.7 


Foregrounds 


7.47 


1.82 


0.948 


Beam 


0.912 


2.00 


4.63 


Point sources 


0.933 


3.75 


5.85 


Final SZ Power 


58.6±17.3 


44.4±50.5 


36.0±50.3 



Table 1. The table derives the SZ angular power spectrum (.{^ + l)Ct/2Ti at the RJ 
end of the frequency spectrum. All values are tabulated in units of /iK^, except for the 
unitless weights. "Raw" and "Final" SZ refer to the spectrum pre and post subtraction 
of the residuals. The Raw SZ spectrum is measured from data, whereas wc simulate the 
contribution from the CMB, instrumental noise, point sources, foreground subtraction 
uncertainties, and beam uncertainties. The residual foregrounds row quantifies the 
residual from the sum of the WMAP estimated foregrounds and the deviation from 
their foreground estimate. These results do not include l//-noise, see the text for 
discussion. The uncertainties are given at 95% c.l. 

161 ± 49 (95% c.l.) in the first bin. This derivation, however, is sub-optimal since it 
does not fully minimize frequency dependent quantities such as instrumental noise and 
foregrounds, but only removes common signals like the CMB. Simulations maintain that 
this method leads to a higher "SZ" signal. Nevertheless, this suggests that the SZ-like 
signal we have detected with our optimal method is not associated with a residual CMB 
that we have not accounted for in simulations related to Table [H 

In binning the SZ spectra into a single multipole bin, we obtain l{l + l)C^/27r = 
55.1 ± 15.6 /iK^ (95% c.l.) for the central multipole of £ = 205, which is heavily 
dominated by the first bin due to its comparatively much tighter error bar. This implies 
a component with 3% the value of the primordial CMB in the power spectrum (or 17% 
in aim temperature). We have failed to find a combination of foreground residuals, beam 
errors, calibration uncertainty, and similar instrumental systematics that will add up to 
the needed correction. 
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Figure 3. CMB mask in Galactic coordinates demonstrating the localization choices. 
The four quarters of the map are numbered accordingly, and the fifth cut is defined 
by the transparent blue band along the Ecliptic plane (±20°). The Galactic plane is 
always masked out with KQ75 in our analysis. 



As a test on the reliability of the signal detected in the first multipole bin, we varied 
the foreground model involving the amplitudes of the synchrotron, free-free, and dust, 
but within the uncertainties allowed by WMAP data we failed to find a model that 
produced a null SZ signal in the first bin. However, a null signal can be obtained if 
the WMAP estimated foreground emission is increased by a factor of two. Our error 
budget includes the full uncertainty in the foregrounds to the extent that foregrounds 
have been established with WMAP data j28| . It could be that our signal is a signature 
of an unknown foreground. If that were to be the case, CMB observations that span 
a wider range of frequencies will become necessary, and we motivate a study with the 
ongoing Planck experiment. 

In order to better understand the physical origin of the measured signal, we repeated 
our analysis for a number of sky cuts. In Fig. [3] we show the CMB mask in Galactic 
coordinates demonstrating the localization choices. Table |2] determines the measured 
signal to primarily lie in the southern hemisphere, at 84.9 ± 24.2 /iK^ (95% c.l.). The 
statistical significance of the signal is equivalent in the two quarters of the southern sky, 
whereas the north-western quarter is completely devoid of an SZ signal. The excess of 
power in the southern hemisphere (Ecliptically aligned) was first detected by Eriksen 
et al (2003) 



36 



and hypothesized to be due to thermal SZ for £ < 10 by Abramo et 
al (2006) |37j. Degree scale anomalies in the CMB have previously been found along 
the Ecliptic plane in Diego et al (2010) 38|, whereas we find the signal in that region 
present at a 2.4(T-level. 
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Large-Scale Signal Localization 

Bin 1 = 10 - 50) Final SZ Power 

Rl (/sky = 0.177) 14.9±32.4 

R2 (/sky = 0.185) 49.2±42.8 

R3 (/^ky = 0.181) 86.8±45.3 

R4 (/sky = 0.173) 58.9±30.8 

R5 (/^j^y = 0.245) 49.2±40.3 

R1+R2 (/sky = 0.361) 34.6±24.8 

R3+R4 (/sky = 0.355) 84.9±24.2 

R1+R2+R3+R4 (/sky = 0.716) 58.6±17.3 

R1+R2+R3+R4-R5 (/sky = 0.471) 57.6±19.6 



Table 2. Rl denotes the north-western quarter of Fig. [3] in Galactic coordinates, R2 
denotes the north-east. R3 the south-west, R4 the south-east, and R5 the Echptic 
plane in a ±20° band. The uncertainties are given at 95% c.l. 



4. Discussion 

If the excess signal we have detected at ten degrees angular scales is indeed SZ, it is 
clear that the signal does not originate from galaxy clusters at high redshift. Such a 
signal would lead to a power spectrum that increases with multipole, rendering an SZ 
signal in the smaller scale bins. A large SZ signal associated with the Extragalactic 
sky is also ruled out by existing arcminute-scale CMB experiments. We do not see the 
signal coming from unresolved point sources or Zodiacal light emission as they do not 



amount to more than a few /iK in strength 38 



One possibility is that the signal is associated with sub-keV gas in the Galactic 
halo, explaining why the signal is primarily at small multipoles. We can try to address 
if such a signal is possible based on requirements on the number density of electrons 
and the existing spectral distortion limits from FIRAS. Since the SZ effect is —2y at 
RJ frequencies, where y is the Compton y-parameter, we can convert the arnplitude 



of the fluctuation power, for a Galactic-like signal with power-law ^ A^~^ [39[] and 
assuming a dispersion in y equal to its average value, to obtain y = (5.1 ± 0.8) x 10~^ 
(95% c.l.). This is consistent with the 95%-level FIRAS constraint of \y\ < 1.5 x 10~^ 
from the CORE satellite 



40 



The electron column density required to render this y-parameter in our Milky Way 

2 

halo can be obtained from Np ^ , where nip is the electron mass, Tg is the 

electron temperature, with other constants having the usual definition. For a 0.3 keV 
electron gas temperature, we find = 1.3 ± 0.2 x 10^2 cm^^ (^^^^ (95% c.l.). This 
electron density is merely approximate due to the assumptions made in deriving y, and 
in presuming all of the measured signal is actually rooted in the SZ. 

The column density of the free electrons is consistent with that suggested in Peiris 
& Smith (2010) [ill to resolve the CMB isotropy anomalies via the kinetic SZ effect. In 
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reality, both the thermal and kinetic SZ effects may contribute to large-scale anomalies 
of the CMB. However, the required column density is in strong tension with pulsar 



dispersion measurements (Taylor & Cordes 1993 |42|]), OVII absorption studies j43 
and the soft X-ray background 4J], which place A'^e ~ 10^° — lO^^cm"^. Given the 
WMAP subtracted foregrounds would need to have been underestimated by a factor of 
two to explain our low-£ signal, this discrepancy suggests that the excess signal may 
originate from further instrumental effects or an unknown foreground. Rather than to 
simply explain away the signal, as we did not find either a simple statistical argument 
or systematic effect, we suggest that further studies are warranted. 
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